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Abstract 



The paper deals with the properties of thermally excited solitons of the 
isotropic spin-S" ferromagnetic chain with nearest-neighbor logarithmic in- 

o 

ly-j ■ teractions. The exact statistical mechanics of the interacting soliton gas is 

developed for the general case (arbitrary S, temperature and magnetic field). 

a ; 

S At low temperatures the model's thermodynamics coincides with that of the 

Heisenberg model. We present analytical approximations of the leading-order 

o ' 

O ■ asymptotic behavior of the energy in three limiting cases: (a) zero field, low 



temperature, classical limit; (b) zero field, T — * 0, S finite (quantum limit); 



(c) zero field, high temperature, classical limit. Cases (a) and (c) are examples 
of a dense gas of [non-topological] solitons; results are in agreement with those 
obtained by the transfer integral method. Case (b) illustrates the behavior 
of a dilute, yet strongly interacting soliton gas; results for the thermodynam- 
ics are very close to (but not identical with) spin- wave and/or Bethe-j4 nsatz 
predictions. 
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I. INTRODUCTION 



Theoretical and experimental work with one dimensional magnets has provided ample 
evidence for the existence and importance of solitons The development of solitonic con- 
cepts and techniques has been characterized by a remarkable degree of cross-fertilization 
between statistical physics, field theory and applied mathematics. Indeed, much of our 
present understanding of the thermal behavior of magnetic solitons results from Mikeska's 
seminal exploitation of the formal equivalence between the anisotropic Heisenberg model 
(with an external field) and the Sine-Gordon continuum. There has been no lack of "direct" 
candidates from magnetism. The isotropic Heisenberg ferromagnet (IHF) has been known 
to be integrable in the continuum limit HQ], and a variant, the Ishimori-Haldane-Fadeev 
ferromagnet (IHFF) is completely integrable on a lattice In both models, the prop- 

erties of solitons at thermal equilibrium, although of fundamental interest, are difficult to 
extract. The reason for this is twofold: (i) the presence of non-topological solitons, whose 
amplitudes can become arbitrarily small and their numbers presumably arbitrarily large, (ii) 
the fact that those non-topological solitons, which perform breather-like internal oscillations 
in addition to their translational motion, actually reduce to linear spin waves in the small 
amplitude limit. The first problem is generic to all systems which carry non-topological 
solitons (e.g. the Toda lattice ||); it means that one must incorporate interactions between 
solitons in order to achieve a valid description of their statistics at any finite temperature 
||. The second problem has been dealt with in the case of Sine-Gordon breathers. There, it 
has been demonstrated (in leading order asymptotics 0]) that if one allows the presence of 
both breathers and phonons, enough of either will be displaced to maintain the total number 
of degrees of freedom. Double counting can in fact be avoided less tediously a priori: exact 
formulations of statistical mechanics in terms of either nonlinear phonons |lTJ or breathers 
12|| yield identical results. 

The present work deals with the thermodynamics of IHFF solitons. At sufficiently low 
temperatures, the model's properties coincide with those of the IHF. Since the [integrable] 
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IHF continuum can be considered to be a good approximation to the IHF lattice under 
the same condition, i.e., at low temperatures, our results in effect provide a description of 
soliton statistical mechanics of the IHF as well. This equivalence has been exploited in a 
recent Letter |13| , where many of our key results have appeared without proof [14| . 



The aims of this paper are: (i) to give a reasonably complete account of the theory of 
[exact] soliton statistical mechanics in the only magnetic system where this has been possible 
to a full extent and, (ii) to support the conclusions of the theory by presenting its application 
to selected analytically tractable limiting cases. Of particular interest in the latter context 



is the treatment of the high-temperature regime (not reported in Ref. [[13]]), which validates 
the concept of a strongly interacting, dense gas of extremely localized solitons. 

The paper is organized as follows: Section II introduces the model and its dynamics. It 
reviews the basic properties of one- and two-soliton solutions and discusses the conservation 
laws which control the symmetries of soliton/soliton phase shifts and are thus central to the 
statistical mechanics. Section III formulates the statistical mechanics in terms of occupation 
probabilities for microstates. Its main result is formulated in terms of a two-dimensional 
integral equation. It will be shown that the exact thermodynamic quantities (soliton density, 
energy, magnetization) are controlled by the asymptotic behavior of soliton quasienergies in 
phase space. Section IV deals with the limiting case of zero-field, high and low-temperature 
asymptotic behavior. In the latter case, it will be possible to derive results for both the 
quantum and the classical regime. Concluding remarks are made in Section V. The rela- 
tionship between thermodynamic properties and integral-equation asymptotics is derived in 
the Appendix. 

II. DYNAMICS 
A. The model 

The IHFF model Hamiltonian 
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H = -2 J5 2 £ ln(i + - gBUSn ~ S) , (2.1) 

describes a chain of iV atoms, each of which carries a spin of length S, and is subjected 
to an external field B along the z-axis. The particular form of the interaction guarantees 
complete integrability |j|-|7| at the classical limit, i.e., as S — > oo, h — > 0, hS — > S'c, J — > 
0, JS* 2 j,g — > 0,g£ — > g c |~5|]. It will, however, be necessary to keep the value of 
as a parameter in the problem, both in order to deal with semiclassical approximations 
of quantum chains and because the equation which describes statistical mechanics at the 
classical limit is singular. 

The classical dynamics of ( |2.1| ) is described by 

^t Sn = W n x Sn (2 - 2) 

where time is measured in units of h/(JS) = S c /j. The fact that the dynamical system 
( p.2| ) is completely integrable has been demonstrated in a variety of ways. Ref. || makes 
use of the gauge equivalence of the dynamics arising from ( |2.2| ) (for B = 0) to the Ablowitz- 
Ladik(AL) |]TB[ version of the discrete nonlinear Schrodinger equation. Refs. ^fj] conjecture 
the existence of such a model by considering a sequence of finite-S" models which are Bethe- 
Ansatz solvable [fL?]]. The complete inverse scattering theory has also been presented |18 
We will follow H and make direct reference to [0 when appropriate. 



B. Properties of single solitons 

The problem of the magnetic field is disposed of first. The physical effect of the field is 
to make all spins perform a uniform precessional motion with angular frequency h = g c B/j. 
This motion can be eliminated by introducing a rotating coordinate frame. Now, since 
the dynamics of the zero-field IHFF ( |2.2| ) can be mapped to the the dynamics of the AL 
model, it is possible to describe any general multi-soliton solution of fl2.2|) in terms of the 
corresponding AL solution, and superimpose the uniform precession due to the magnetic 
field at the end. This simplification is of particular importance for the statistical mechanics, 
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since it means that all relevant properties of interactions between solitons are controlled 
by the values of their dynamical parameters for vanishing magnetic field; it does not imply 
however that thermodynamic properties have a trivial dependence on the magnetic field (c/. 
Section III). 

The form of the single soliton solution (cf. demonstrates that all physical soliton 
properties can be described in terms of two parameters, w and k. Thus, the soliton has a 
spatial extent (in units of the lattice constant £) equal to 1/ {2w), it moves with a translational 
velocity U(w,k) = sinh2wsin2fc/(2w), and performs an internal oscillation of frequency 
Q(w, k) = 2kU(w,k) + 2(cosh2wcos2A; — 1) + h; furthermore, one can evaluate the two 
obvious integrals of motion, namely, the total magnetization 

sinh 2w 

M(w, k) = -2 — -, 2.3 

v ' ; cosh 2w- cos 2fc' V ' 

(in units of S c ), and energy 

E{w,k) =8w - hM(w,k) (2.4) 



(in units of j) carried by a single soliton. The soliton canonical momentum P, defined [19 
via {dP/dE) M = l/U, (in units of S c /£), determines the correct measure to be used in 
statistical mechanics: 

dPdM = 8 [ P ' M } dwdk = f- ( dwdk = 32- — dwdk (2.5) 

d{w,k) U\dkj w (cosh2w-cos2A;) 2 v ' 

In the limit w — > 0, the properties of a soliton characterized by (w, k) reduce to those of a 
plane spin wave with wavevector q = 2k and frequency 4sin 2 (g/2). 



C. Multisoliton solutions: soliton-soliton phase shifts 

The simplest way to evaluate the asymptotic space shifts resulting from the interaction 
of two solitons is to go back to the AL formulation of the inverse-scattering transform, write 
down the two-soliton solution, and examine its asymptotic form as t — > ±oo. The calculation 
is tedious but straightforward. Here, we only quote the result and make a few comments. 
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The phase shift which a soliton (w, k) experiences due to its collision with a soliton (w', k') 
is given by 



A(w,k;w',k') = — In 
2w 



It follows that the sum rule 



cosh 2 (w' + w) — cos 2 (A; — k') 
cosh2(w' — w) — cos2(A; — k') 



(2.6) 



2ix 

dk A(w, k; w', k') = — mm(w, w') (2.7) 

w 

is always satisfied. Furthermore, the soliton/soliton phase shifts obey the following symme- 
try (generalizable to any multisoliton solution characterized by {w\, k\}): 

]T w x A(w x ,k x ;wy,k y ) =0 (2.8) 

A^A' 



D. Phase shift symmetries and conservation laws 

Symmetries of the type ( |2.8| ) are generic to soliton-bearing systems. They reflect the 
existence of a low-order integral of motion / = J2x w\U\ , (where U\ is the velocity of the 
Xth soliton) and imply that the [conserved] quantity J2x w \ travels at a constant speed. 
Inspection of the formula for the group velocity shows that the required conserved quantity 
is AL's Ci, or, in the spin language, 

I = y^Xn = ^ n ' X - (2.9) 

n n (1 + S n ■ 5 n+ i)(l + S n ■ S n —i) 

and reflects the invariance of the original spin Hamiltonian under infinitesimal rotations of 
each spin S n around an axis defined by the spin-space gradient 

V S,/ = V 5„(Xn+l + Xn + Xn-l) ■ (2.10) 

Note that generally {i.e. for h 7^ 0) w is proportional to the exchange part of soliton energy. 
It is possible to show that the center of exchange energy of the entire IHFF lattice travels at 
a constant speed, independently of any soliton considerations. Given an infinite chain with 
any type of isotropic nearest neighbor spin interaction f(S n ■ S n+ \), consider the quantity 
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X = -Y,ng(S n -S n+1 ), (2.11) 

L n 

where g is such that the denominator Z = J2n[g(S n ■ S n+ i)} is a constant of the motion. The 
derivative dX/dt can be written down by using ( |2.2| ); it turns out that, in order for it to 
remain constant under the Hamiltonian flow, (i) g must be proportional to / (and thus Z 
proportional to the exchange part of the Hamiltonian) and, (ii) the quantity 

^2g'(S n ■ S n+ i)g'(S n ■ S n -i)S n ■ (S n+1 x 5 n _i) (2.12) 

n 

must be a constant of the motion, as indeed is the case with the logarithmic (but not with 
the Heisenberg) interaction. Thus, the existence of the constant of motion / guarantees a 
very special status to the exchange energy of the IHFF model (analogous to that of the mass 
in a lattice model ||). 



III. SOLITON THERMODYNAMICS: EXACT RESULTS 



The general part of the formalism used has been presented in [20|, in the context of 



the Sine-Gordon equation; the reader is referred there for details. The physics involve the 
minimization of the classical free energy functional with respect to the occupation probability 
of each state. The procedure is self-consistent since the density of states is itself a functional 
of the occupation probabilities, due to the pairwise additive interactions ( [2.6|) which restrict 
available phase space. The thermodynamics derived within this framework is in principle 
exact, with a single caveat: because solitons with a very large characteristic length are in 
practically indistinguishable from linear excitations (spin waves, cf. above) one expects 
that a correct statistical mechanical treatment will ultimately eliminate one or the other. 
The situation is analogous with the famous case of the Sine-Gordon breather. It appears 
therefore expedient to formulate the theory in terms of solitons alone (cf. ||21||). 

The thermodynamic potential (free energy) of the interacting soliton gas is given by 

-j3F = J dTR (T)n(T) (3.1) 



where (3 = JS 2 /ksT is the dimensionless inverse temperature, V is shorthand for the point 
(w, k) in the soliton parameter phase space, and the phase space element dT is defined by the 
measure ( |2.5| ). The density of states -Ro(r) refers to the non-interacting soliton gas (with 
our choice of units (cf. Section lib) equal to N/(2ti)S 2 ) and the occupation probability 
n(r) = exp(— (3e) defines a characteristic [field- and temperature-dependent] quasiparticle 
energy e(T). The latter is determined by the integral equation 

/3e(r) = pE(r) + - 1 - I dr'ACr', r) e - /3e(r ' ) . (3.2) 

Z7T J 

A further quantity of interest is the density of available states at any given point in phase 
space expressed as fraction of Ro(T), i.e., R(T) = R(T)/Rq(T). In the IHFF case, equations 
J) - Q take the form 

_Pf=±^ r dk r dw - e"^.*) (3.3) 

7r Jo Jo (cosh 2w — cos 2k) 2 



and 



Q /*7T /*00 

Pe{w, k) = (3E(w, k) + -S 2 / dk / dw' 

7T JO JO 



In 



(cosh 2ui'— cos 2k') 2 

cosh 2(ui'+ui)— cos 2(k—k') 
cosh 2(ui'— ui)— cos 2(k— k') 



e-P< w '> k '\ (3.4) 



respectively, where / is now the free energy per site. Furthermore, 



E iWlk )='( d ^ k) ) (3-5) 



0h 



The derivation is straightforward, and involves use of (|2.5| ) and ( |2.6|) . Some special attention 
must be paid to the analogy of ( |3.5|) with Eq.(19) of p0[] ; in this case, the crucial symmetry 
refers to the exchange part of the energy of a free soliton (cf. above); hence the denominator 
E (w) = 8w in (p.5|). Moreover, the derivative with respect to (3 must be taken at constant 
Ph. 

The thermodynamic functions we are interested in, i.e. energy, soliton density and mag- 
netization can be obtained by differentiation of ( ^.3[ ) with respect to (3, fi and h respectively. 
The soliton chemical potential must vanish, since the soliton density is not a priori fixed. 
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The details are presented in the Appendix. The important feature is that we do not require 
detailed knowledge of the solution of fl3.4| ) at all points in the two-dimensional region defined 
byO<u><oo,0<A;<7r. The thermodynamic behavior of the soliton gas is controlled by 
the asymptotic properties of ( \3.$ ) in the limits w — > and w —>■ oo. In particular, the free 
energy is given by 



Ztt Jo 



and the soliton density per site by 



where 



n s = — / dk 



' dV 



I3h 



V{w,k) 



d(3e(w, k) 
dw 



(3.6) 



(3.7) 



(3.8a) 



= lim V(w,k) = 8/3, 



(3.8b) 



Vq = lim V(w, k), and 



(3.8c) 



Coo(fc) = lim ((3e(w, k) — (3E(w, k)) 



(3.8d) 



It follows that the internal energy u and magnetization m per site can be obtained by direct 
differentiation of (|3.6| ) with respect to (3 and h respectively, i.e., 



u 



^d(3j h 2nJ 



dk 



'ocuky 



(3.9) 



m 



'w\ ._j_r dk 

dh L 2tt Jo 



dh 



(3.10) 
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IV. SOLITON THERMODYNAMICS: ANALYTIC APPROXIMATIONS 



A. The classical, zero-field, low temperature case 

At low temperatures we expect thermodynamic properties to be dominated by solitons 
whose exchange energies are lower than 1/(3, i.e., with parameter values w <C 1. If, in 
addition, k > w, it is possible to approximate the exact phase shift Q2.6| ) by 

A(w, k;w' , k') = — min(w, w') S(k — k') . (4.1) 

w 



Note that the above approximation preserves the symmetry (|2.8|) of the soliton/soliton phase 



shift and satisfies the sum rule ( J2.7|) , so we can reasonably expect it to exhibit at least the 
salient features of the exact solution of (|3.4j) . Introducing (|4.1|) in Q3.4|) , it is possible to 



decouple the statistics of solitons with different A;— values, and transform the two-dimensional 
integral equation ( |3.4|) to a set of one-dimensional integral equations in which k appears as 
a parameter: 

Pe(w, k) = (3E(w, k) + 8S 2 f°° ^ -W,.) (42) 

Jo (w' z + sin ky 



The integral equation ([4.2|) can be transformed to a second-order differential equation for 

x = (3e{w, k): 

( p. x e -x-f3hM{t,k) 

dP = " 8 (1 + at*)* ' (43) 

where a = sin 2 k/S 2 and we have used a rescaled independent variable t = SwJ sin 2 k, in 
terms of which M(t, k) = —2t/ (1 + at 2 ) / S. Eq. (|4.3|) must be solved subject to mixed initial 
and final conditions, i.e., x(0) = 0,dx/dt(oo) = f oo = 8((3/S) sin 2 k. (Note the rescaling 
with respect to of ( p.8£|) ) 

The remainder of this Section specializes to the case h = 0. 

In the classical limit, S — > oo,/3 finite, the quantities 

Voo = 8(3 sin 2 k/S (4.4) 
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and a become infinitesimally small. Eq. ( f4.3|) can be further simplified to 

— -8 ~ x 
dt* ~ ' 

which, for the given boundary conditions, has the exact solution 
,(^ 2 l n {±sinh(^ + si n h-(^))}. 



(4.5) 



(4.6) 



Two remarks are in order here. First, that exactly the same equation (|4.5| ) controls the 
classical limit of the Sine-Gordon breather gas Second, that if one regards ( |4.3|) , with 
h = 0, as a mechanical system, (|4.5|) represents its "autonomous" limit. The asymptotic 



properties of (^4.6|) are particularly simple (c/. definitions above): 



v = y/lQ + v, 



2 

oo i 



(4.7) 



and 



Coo = 21n 



VQ + Vc 



(4.8) 



It follows that 



(4.9) 



In the classical limit under consideration, This means (cf. Eq. (577)) that the magnetic 
lattice will be equally populated with solitons of all k, producing a total density of 1/2 
solitons per site, as expected from counting the classical degrees of freedom. The internal 
energy density can be found by inserting ( |4.9|) and (|4.4j) in (|3.9| ) to be equal to k B T, as 
demanded by the equipartition theorem. 



B. The finite S, zero- field, f3 — > oo limit 

The core of the approximation scheme developed in the preceding subsection remains 
valid away from the classical limit if the temperatures become extremely low. More precisely, 
for S finite and (3 large, consider the range of k defined by 
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8(13 /S) sin 2 k < A 



(4.10) 



where A ^> 1 but finite. As (3 — > oo, this range is bounded by a characteristic k c oc l/yf/3. 
For such fc's, the parameter a is smaller than In the limit /3 — > oo this allows us to 

substitute 



by (4-5) . On the other hand, values of k outside this range are (in the limit 
(3 — > oo) irrelevant. The argument for this is as follows: such A;- values lead to > >1> 1; 
numerically large values of lead in turn to comparably large values of Vq, in the context 
of both ( |4.5| ) and ( f4.3| ). In the former case this can be proved by inspecting ( |4.7| ). Moreover, 
it is in both cases clear that a ^> 1 can only come from a «o > 1. This however means 
that x(t) grows very rapidly from its initial value x(0) = 0. As this happens, the right 
hand sides of both (|4.5|) and (|4.3j ) vanish, dominated by the exponential factor. Since the 
"acceleration" vanishes before it has any time to act, the "velocity" stays at its initial value, 



I.e., « Vq, 



dvo 

dVr^ 







(4.11) 



and therefore such solitons are not thermally excited (c/. fl3.7|) ). As a consequence, in the 
limit (3 — > oo, it is legitimate to use (4.5) in lieu of (|4.3|) for all k, and thus exploit the 
properties flE7D-(pESD of the solution flO|) . 

Inserting ( |4.10| ) in ( p.7|) , and approximating sink by its argument, we obtain the leading- 
order contribution to the soliton density, 



n, = —\ — - / dx I 1 — 



1/2 



7i V 2/3 Jo 



X 



X 1 



JS 



(4.12) 



The corresponding result for the energy density, obtained by making use of ( [4.10 ) in ( (4.9|) 
is 



u = 2- 



3/2 



(4.13) 



A few remarks are in order here. First, Eq. ( |4.12[ ) means that the soliton gas is dilute. 
Not all solitons are excited. This is a characteristic of the quantum regime. In the classical 



12 



limit (c/. previous and following subsection) solitons must exhaust all degrees of freedom. 
(It should be remembered that no provision for linear excitations has been made.) Second, 
solitons carry an average energy of IksT each, as suggested by equipartition. Numerical 
solution of ( |4.3j ) indicates |H| that this result is only valid in the limit (3 — > oo. Third, the 
reason the approximation ( |4.1[ ) works is that the conditions for its validity are met by thermal 
solitons; characteristic w- values of are much lower than the typical k values of 0(k c ) 
(quantum case), or 0{\) (classical case, preceding subsection). Fourth, the numerical value 
of the prefactor of the energy, 0.3814, is close to the spin-wave result of 0.3684. Comparable 
results have been obtained in the 5=1/2 case in terms of the Bethe-Ansatz |2^j2^]. In 
view of the soliton-magnon duality at low amplitudes, this is not unexpected. 



C. The classical, zero- field, high temperature limit 

In the high temperature regime, (5 < 1, narrow solitons which carry large amounts of ex- 
change energy, i.e., with values of w > 1 become predominant. As long as the characteristic 
parameters w,w' of two such solitons do not come very close to each other, the hyperbolic 
cosines in the exact soliton/soliton phase shift ( |2.6p can be approximated by exponentials, 
leading to 

2 

A(w, k; w', k') = — mm(w, w') . (4.14) 
w 

The above approximation is the exact opposite of ( |4.1|) : instead of a 5 function in fc-space, 
which decouples soliton modes corresponding to different fc's, (|4.14| ) is independent of k, k', 



ensuring that all soliton modes couple to each other. It should be further noted that (|4.14|) 
satisfies the sum rule (|2.7|) , just like the low-T approxination ( |4.1| ). 
Substitution of (|4lD in Q Y ields 

Pe(w, k) = 8f3w + / dk' / dw'e- 4w ' min(w, w')e-^ w '' k,) , (4.15) 

where again the hyperbolic cosine in the denominator of ( ^.4[ ) has been approximated by an 
exponential, to ensure consistency with ( (4.14[) . It should be emphasized that these approx- 



13 



imations are crude, and that they cannot be expected to yield accurate results for values of 
w <C 1. The point is that such low energy solitons are irrelevant at high temperatures. 

In accordance with this argument, we temporarily exclude soliton modes with very low 
values of w from any further consideration by introducing a lower cutoff in ( [4.15 ). This will 



be useful in controlling singularities, although the value of the cutoff will be set to zero at 
the end of this calculation. 

Since the right-hand side of ( |4.15|) is /c-independent, e can only depend on w; Integration 



over k' can be performed trivially, and reduction to an o.d.e. follows the steps of the previous 
subsection. The result is 

4^ = -U8S 2 e^ x+4 ^ . (4.16) 

dw l 



The most general solution of Eq. ( ^4 . 1 6| ) consistent with the condition lim^^oo dx/dw = 
Voo = 8(3 is 

x(w) = 21n|^e- 2M, 5sinh(2Aw + 6)|, (4.17) 

where A = 1 + 2(3. We do not know the value of b, since we have no knowledge of the exact 
behavior of x(w)as w — > 0. The following procedure will be adopted instead: b will be set 
equal to zero, and the value of the cutoff will be fixed by demanding x(wq) = 0. This leads 
to Wo = 1/(8S), which vanishes at the classical limit. The solution (|4.17|) has the asymptotic 
properties 

VfW(^) = -4+16S (4.18) 



and 



Coo = -2 In f^-^] +2 In 5. (4.19) 



Both (|4.18|) and ( |4.19 ) are valid for all k. In the limit S — > oo, V diverges, but 
lim^oo dVQ /dVoo = 0, and thus, from ( ^.7[ ) the soliton density must equal 1/2. Fur- 
thermore, from ( |3.9|) and fl4.19|) , the internal energy per site 
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can be calculated, in agreement with the exact result obtained by the transfer integral 
method [03 . 



Although this completes the derivation of soliton thermodynamics at the classical high 
temperature limit, it is instructive to examine the basic features of the analytical solution 
( 4.17 ). For example, using (|4.17|) in the alternative definition of the soliton density (c/. Ref. 
2(1), allows us to examine the u>-distribution of thermal solitons (with the understanding 



that no conclusions should be drawn about details in the region 0<w<l): 

poo i ( 2\w 1 

n s = 2X dw <^ , - r -1} . (4.21) 

Jo sinh 2 (2Aw) \tanh(2Aw) J v ; 

Eq. ( |4.21| ) implies that solitons which are thermally excited have values of w ~ 1/A. In 

other words, typical thermal solitons cannot become narrower than a lattice constant. Since 

each such soliton carries an energy of order j, the total energy per site is of order unity at 

high temperatures, as expressed by (|4.20|) . 



V. CONCLUDING REMARKS 

We have presented an "advanced phenomenology" of the interacting soliton gas for an 
isotropic ferromagnetic chain. The theory is in principle exact but requires solving a two- 
dimensional integral equation, i. e. in effect has the same degree of complexity as the quan- 
tum mechanical Bethe-Ansatz. In the case of vanishing magnetic field it has been possible 
to develop valid approximations in the high and low temperature regime and derive an- 
alytic expressions for thermodynamic quantities. The results presented provide detailed 
information about the manner in which nonlinear modes accommodate themselves at finite 
temperatures. 

In the low temperature regime, the fortunate occurrence of a gapless spectrum allows a 
semiquantitative understanding of the quantum, T — > 0, limit of the Heisenberg chain with 
essentially classical techniques. The behavior of the soliton density suggests the presence 
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of a "dilute" soliton gas. The quotation marks imply that the term should be used with 
some caution. It refers to typical distance between the centers of two solitons. At low 
temperatures, this is indeed large ( of order T -1 ' 2 ); however, the spatial extent of a thermal 
soliton is of order 1/T, i.e., much larger. There is therefore significant overlap between 
solitons, suggesting the presence of strong interactions. It is interesting to digress at this 
point and consider the analogous situation in the Toda lattice. There, the typical intersoliton 
distance at low temperature is of the same order as the spatial extent of thermal solitons, 
i.e., T- 1 /' 3 . 

It appears that a properly formulated soliton phenomenology has a range of validity 
which far exceeds the limits imposed by the original conceptual framework p5| , |26[| . In the 



particular case of the quantum low-temperature regime, marginal uncertainties in the Bethe- 
Ansatz make it presently impossible to decide whether the soliton result is asymptotically 
exact. It is however clear, that soliton theory presents a better alternative to spin waves: 
the magnetization obtained from the soliton theory vanishes at any finite temperatures (c/. 
the divergent magnetization obtained from spin- wave theory). 

At high temperatures it appears, that in spite of the crudeness of our approximation, 
the theory again describes leading order corrections. The exact agreement of ( |4.2(J| ) with 



the transfer integral result should be regarded as fortuitous. In this case the soliton gas is 
dense in the traditional sense. Intersoliton distances and soliton spatial extent are both in 
the range of the lattice constant. 

APPENDIX: 

In order to derive ( |3.6|) , we subtract the "bare" part of the soliton energy, /3E(w, k) from 
both sides of (|3.4[) and integrate both sides over all k. The result is 



/ dk{Pe-(3E(w,k)} = 8S 2 / dk dw'^—— f — — 2 min( W , w')e^< w > k > , 

Jo Jo Jo (cosh 2u>' — cos 2k' ) z 

(Al) 
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where we have made use of the sum rule (|2.7| ). In the limit w — > oo, the integrand of the 
left-hand side of ( |AI| ) is C^, whereas the right-hand side coincides, except for a prefactor 
l/(27r), with -Pf (cf. Eq. U). This proves fl3j). 

In order to derive (|3.7|), we differentiate both sides of ( |3.4j ) with respect to w. In the 
limit w — > 0, this yields 



A/3h 8S 2 r* „ /"°° , , 1 4sinh2w' 



= + - ^ o; + — Tdk Hdw' 

1 — cos 2/c 7T jo jo 



-/3e(V,fc') 



1 — cos 2/c 7T jo jo (cosh 2k/ — cos 2k') 2 cosh 2k/ — cos 2(k' — fc) 

(A2) 

Differentiating both sides of ( |A2|) with respect to at constant (3h, we obtain 

(dv„(k)\ _ 1 _^_f Jk i- dw , W — - 4sinh2m ' - it {w ', vyr****) , 



dVoo J /3h it Jo Jo (cosh 2k/ — cos 2//) 2 cosh2K7' — cos2(// — k) 

(A3) 



where we have used ( |3.5| ) for the density of states R. Eq. can now be obtained by 
integrating both sides of (|A3[) over k. 

A further by-product of ( |A3| ) is that, in the limit k — > 0, the integrand in the right-hand 
side expresses the total magnetization in soliton component form. This implies that the 
total magnetization per site (in units of S c ) can be alternatively expressed in terms of the 
asymptotic properties of (|3.4j) , as 

m = hm " v ; . A4 

The above alternative can be useful in situations where the integral ( |3.10| ) cannot be calcu- 
lated exactly. In all limiting situations which we have treated analytically (cf. Section IV), 
^) guarantees a vanishing magnetization. 
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